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1.  Introduction 

The  application  of  the  Wald  SPRT  to  life  testing  of  equipment, 
having  an  exponential  distribution  of  the  interfailure  times,  has  been  in 
practice  from  the  fifties ZacRs- {1971-?  pp .  46b-]->.  ~*The  problem 
discussed  in  the  present  paper  is  that  of  the  exact  computation  of  the 
distributions  of  the  stopping  times,  and  the  determination  of  confidence 

"Til  f’  1 

limits  for  the  meantime  between  failures  (MTBF),  6,,  after  stopping.  -  'i* 

More  specifically,  we  consider  a  life  testing,  or  reliability  acceptance 
test,  in  which  n  identical  systems  are  put  on  test.  VJhenever  a  failure 
occurs  it  is  instantaneously  repaired,  the  time  of  failure  is  recorded, 
and  the  system  continues  to  be  tested.  The  total  number  of  failures 
of  the  n  systems,  in  the  time  interval  (0,t],  is  a  Poisson  process, 

X^(t),  with  mean  nt/6,  where  6  is  the  MTBF.  In  Section  2  we  present 
a  modified  Wald  SPRT,  for  the  determination  of  the  stopping  time,  and 
the  decision  of  whether  to  accept  or  reject  the  hjrpotheses 
Hq:  6  >  6q  vs.  Hj^;  6  <  6^^,  for  some  0  <  6^^  <  6^  <  In  Section  3 


we  define  the  distributions  of  the  stopping  times  at  acceptance  or  at 
rejection,  and  the  distribution  of  the  actual  stopping',  time.  In  Section 
A  we  show  how  one  can  determine  confidence  intervals  for  6  after  stopping. 
The  formulae  on  which  the  algorithm  for  the  determination  of  the  distribu¬ 
tion  of  the  stopping  times,  the  acceptance  probabilities  and  average  test 
length  are  given  in  Sections  5  and  6.  A  BASIC  program  for  the  computa¬ 
tions  is  given  in  the  Appendix. 

The  problem  of  determining  confidence  intervals  after  sequential 
stopping  is  discussed  by  Slegmund  (1978,  1985)  and  Wijsman  (1981). 

Siegmund  shows  how  one  can  apply  results  from  renewal  theory  to  obtain 
approximations  to  the  required  distributions.  By  utilizing  the  specific 
properties  of  the  Poisson  process,  as  a  non-decreasing  process,  of 
unit  jumps,  we  have  developed  recursive  equations  for  exact  computations. 
The  algorithm  is  practical  with  the  available  modern  equipment.  It  can 
be  easily  executed  with  a-  desk  top  computer. 


Boundaries  for  Reliability  Testing 


In  Figure  1  we  present  the  stopping  boundaries  of  a  modified 

SPRT.  The  modification  of  the  SPRT  is  in  the  introduction  of  an  upper 

bound,  k  ,  for  the  number  of  allowed  failures.  We  assume  that 
s 

k  >  krt  +  k, .  Let  X  (t)  denote  the  total  number  of  observed  failures 
s  -  0  1  n 

of  n  units  on  test,  in  the  time  interval  (0,r].  The  rule  is  that  the 
experiment  is  terminated  when  the  process,  X^(t),  either  crosses  the 
upper  (rejection)  boundary,  or  intersects  the  lower  (acceptance) 
boundary.  The  upper  (rejection)  boundary  is  given  by 


!k^  +  bt,  if  0  <  t  <  t 
s  _  _ 

where  k^  is  a  positive  (intt  er)  intercept,  and  t  =  (k^  -  k^)/b.  Tl. 


slope  b  is  determined  to  achieve  certain  characteristics,  which  will  be 
discussed  in  the  next  section.  The  lower  boundary  for  acceptance  is 
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given  by  the  line 


+  bt,  0  <  t  <  t  ;  (2.2) 

an  integer.  The  value  of  t,  =  k„/b 

Kq  U 

ceptance.  By  increasing  the  value  of 
the  minimal  acceptance  time  and  thus 
decrease  the  probability  of  accepting  the  product.  By  Increasing  the 
value  of  we  decrease  the  probability  of  rejecting  the  product,  and 
by  increasing  the  value  of  the  slope  b  we  increase  the  probability  of 
accepting  the  product. 

3*  Distribution  of  Stopping  Times  and  the  Characteristics 
of  the  Sequential  Procedure 

We  distinguish  between  two  stopping  variables,  stopping  while 
crossing  the  acceptance  boundary,  and  stopping  while  crossing  the 
rejection  boundary.  We  denote  the  stopping  time  at  acceptance  by  T., 

A 

and  the  other  one  by  The  life  testing  terminates  at  time  epoch  T, 

which  is  the  smaller  of  and  Tj^,  i.e.,  i  =  min{T^,  Tj^).  Obviously, 

either  T  =  or  t  =  Tj^.  Let  (t)  =  f  cumulative 

distribution  (CDF)  of  T^,  when  the  MTBF  is  9.  Let  FQ^(t)  be  the 

corresponding  CDF  of  Tj^,  and  F0(l)  that  of  t.  Since  one  terminates 

either  with  i  or  with 

A  R 

(3.1) 

and  since  T  <  t^,  where  t^  =  t  -  1/b,  ^0(ig)  *=  1  •  On  the  other  hand 

R 

^0  (^s^  ^  ^0  (^s^  Notice  also  that  t  <  T.,  and  t.  can 

Kq  “A  A 

accept  only  the  values  t  .  .  =  t,  +  1/b,  1=0, 1,2 . k  -  1. 

^  *^0  ^  *^0  ® 

Accordingly,  Fg  (t)  =  0  for  rll  t  <  t,  .  The  distribution  of  x  is 

Kq  a 

discrete,  with  jumps  at  t  ,  of  size 

Rq  +  1 

-  4  - 


■  -'*0 

where  t  =  (k  +  k-)/b,  and  k_  Is 
s  J  u 

is  the  minimal  time  required  for  ai 
kg  (with  a  fixed  slope  b)  we  defer 


The  value  of  F  ^(t)  at  t  ,  =  t.  +  (k  -  l)/b  is  the  probability  of 

O  s  *^0  ® 

acceptance.  We  denote  this  probability  by  71(6).  The  distribution  of 
^R’  ^0  other  hand  is  continuous.  Explicit  formulae  for 

Fq  (t)  and  F^  (t)  are  given  in  Section  5.  In  Table  1  we  present  numerical 

values  of  these  distributions,  for  the  case  of  N  =  20,  =  6, 

k^  =  11,  k^  «  50,  b  ■=  20  and  0  =  3,500  [hr].  We  see  in  this  table  that 

the  acceptance  probability  with  the  specified  boundaries  is  tt  (3,500)  =  .A121. 
The  10-th  percentile  of  F„(t)  is  t  ,_  ■=  .  75  x  3,500  =  2,625  [hr]. 

The  median  is  T  =  2.13  x  3,500  =  7,455  [hr]  and  the  90-th  percentile 

is  T  gg  =  2.55  X  3,500  =  8,925  [hr].  The  expected  stopping  cime,  or 

average  test  length  (ATL)  can  be  easily  computed  by  the  approximation 
formula 

ko+kg-1 

ATL(e)  =  i  (i  -  ■|)[Fg(t^)  -  Fg(t^_^)],  (3.3) 

where  =  i/h,  1  •=  ],2,...,k^  +  kg  -  1,  and  tg  ^  0.  A  formula  for  the 

exact  computation  of  ATL(0)  is  given  in  Section  6.  In  Figure  2  we 
present  the  graphs  of  the  acceptance  probability  fui  ction,  tt(6),  and  the 
functions  ATL(6)  and  the  median  stopping  time  T  ^(6). 

We  see  that  the  average  test  length  attains  a  maximum  value  of 
1.9  [3,500  hr]  when  6  =  3,500  [hr].  On  the  other  hand,  when  6  <  2,500 
[hr]  or  6  >  4,500  [hr]  then  ATL(0)  <  3,500  [hr],  which  is  almost  50% 
savings  in  total  time  on  test.  There  is  some  discrepancy  between 
ATL(e)  and  T  ^(6).  When  6  <  2,500  [hrj  or  6  >  3,750  [hr],  ATL(0)  >  T  ^(6). 


•T'  »"T  TTJ’TCT^^*WW>lWPaB"-7'V''-r  ir^TTH  T 


w-^TrT/ir-fffr  TTywv^-’ 
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Table  1 

The  distributions  of  T^,  and  T, 
for  •=  11,  =  6,  k  =  50,  b  ■=  20. 

X  US 

6  ■=  3500,  T  is  in  units  of  3500  hr 


T 

F^(t) 

F^(T) 

Fg(t) 

0.05 

O.OOOo 

0.  10 

0.  OOt'O 

0.  15 

0.0000 

0.20 

0.  Oi.lOO 

0.25 

0. ooOl 

0.  JO 

0-0002 

O. 0025 

0. 0027 

0.35 

0.0004 

0. 0080 

0.0084 

0.40 

O. 0008 

0.0160 

0.0168 

0.45 

0.0013 

0.0260 

0.0273 

0.50 

0. 0020 

0.0373 

0.0394 

0.55 

0.0029 

0.0496 

0.0525 

0.60 

0.0040 

0.0623 

0.0663 

0.65 

0.0053 

0.0753 

0.0806 

0.70 

0.0068 

0.0883 

0.0952 

0.75 

0.0086 

0. 1013 

0. 1099 

0.80 

0.0105 

0. 1141 

0. 1246 

0.85 

0.0126 

0. 1266 

0. 1392 

0.90 

0.0148 

0. 1389 

0. 1537 

0.95 

0.0172 

0. 1509 

0. 1681 

1.00 

0.0198 

0. 1625 

0. 1823 

1.05 

0.0225 

0. 1738 

0. 1963 

1.10 

0.0233 

0. 184/ 

0.2100 

1.15 

0. 0282 

0. 1953 

0.2236 

1.20 

0.0313 

0.2056 

0. 2369 

1.25 

0.0344 

0.2156 

0.2500 

1.30 

0.0376 

0.2252 

0.2628 

1.35 

0.0409 

0.2346 

0.2755 

1.40 

0.0442 

0.2437 

0.2878 

1.45 

0.0476 

0.2524 

0.3000 

1.50 

0.0510 

0.2610 

0.3120 

1.55 

0.0545 

0.  2692 

0.3237 

1.60 

0.0579 

0.2772 

0.3352 

1.65 

0,  0615 

0. 2850 

0.3465 

1.70 

0.0650 

0.2926 

0.3576 

1.75 

0.0686 

0. 2999 

0. 3685 

1.80 

0.0721 

0. 3070 

0.3791 

1-05 

0.0757 

0.3139 

0.3B96 

1.90 

0.0793 

0.3207 

0.3999 

1.95 

0.0820 

0.3272 

0.4100 

2.00 

0.0935 

0.3336 

0. 4271 

2.05 

0.  n  11 

0.3398 

C.4509 

2.  jO 

0. 1355 

0. 3450 

0. 4814 

2.  15 

0.  1063 

0. 3517 

0.5180 

2.20 

0.  2026 

0. 3574 

0.5601 

.>  oc 

0..  243  7 

0. 3630 

0.6068 

2 . 30 

0.2883 

0. 3685 

0.6568 

2.  35 

0. 3350 

0. 3738 

0. 7087 

2.  40 

0.3820 

0. 3790 

0.7610 

2.  45 

0.4277 

0. 3B4C 

0.8118 

2.50 

0.4704 

0.3890 

0.8594 

2.55 

0.5083 

0.3938 

0.9021 

2.60 

0. 5400 

0.3985 

0.9385 

2.65 

0.5643 

0.4031 

0.9674 

2.  70 

0.5804 

0. 4077 

0.9880 

2.75 

0.5879 

0.4121 

1  .  C>00<  1 
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This  reflects  that  the  distributions  F  (t)  are  positively  skewed  in 

y 

those  ranges.  On  th;  other  hand,  if  2,650  <  6  <  3,750  then 

T  ^(9)  >  ATL(9),  corresponding  to  negative  skewness  of  F^Ct).  The 

graph  of  IT  (6)  shows  that  the  probability  of  accepting  a  product  whose 
MTBF  is  less  than  3,000  [hr]  is  not  exceeding  .1.  On  the  other  hand, 
if  the  MTBF  is  greater  than  4,500  [hr]  then  11(6)  >  .9. 

4.  Confidence  Limits  for  the  MTBF 


In  Figure  3  we  present  curves  for  the  determination  of  the 
confidence  limits  for  the  MTBF,  as  a  function  of  the  stopping  time  T. 
If  T  «=  we  use  the  curves  at  the  acceptance  region.  The  curve 

designated  by  LCL  is  actually  the  graph  of  the  20-th  percentile  of 

•  o 

the  distribution  of  T^,  Fg  (t),  as  a  function  of  6.  We  obtain  this 

graph  by  computing  the  distribution  Fg  (t)  at  various  values  of  G 

greater  than  3,175  [hr].  As  seen  in  Figure  2,  71(3,175)  .2.  Thus, 

the  graph  of  2^®)  at  6  =  3,175  is  located  at  »  2.75  [3,500  hr], 

which  is  the  largest  value  that  T  can  assume.  We  find  the  lower 

A 


confidence  limit  for  6,  at  confidence  level  of  Y  ^  -S,  by  reveise 
Interpolation  in  Figure  3.  If  for  eiiample,  in  an  actual  life  testing 
we  terminate  with  7^  =  .75  [3,500  hr],  we  read  in  Figure  3  the  value  of 

0  which  yields  7  (9)  =  .75,  namely  6  =  4,000  [hr].  The  lower 

A  I  «  ^ 

confidence  limit  is  ^  «  4,000  [hr].  If  7  .5  [3,500  hr] 

•  8  A 

then  ^  (.50)  •=  5,000  [hr].  Notice  that  the  corresponding  lower 

confidence  limits  for  Y  •=  .9  are  read  from  the  graph  labeled  LCL  ^ 


to  be  0  g(.75)  =  3,500  [hr]  and  G  ^(5)  =  3,875  [hr].  If  7  “  7 

that  is  the  test  itself  rejects  the  product,  upper  confidence  Intervals 
for  6,  at  confidence  levels  of  y  *  -b,  .9  and  .95  can  be  determined 
by  using  the  curves  at  the  rejection  region  of  Figure  3,  If,  for 


8 


MTBF  [Hr] 


example,  T  ■=  1.25  [3,500  hrj ,  we  read  from  the  graphs  that  the  upper 
confidence  limit  at  Y  *  .8  is  6  g  *  2175  hours.  If  assumes  larger 
values  there  is  an  indication  that  6  is  larger.  For  example,  if  Tj,  * 

2.00  [3,500  hr]  then  the  confidence  interval  is  (2,450,  3,175). 

5.  Computing  the  Distributions  of  Stopping  Times 

Let  X^(t)  designate  the  total  number  of  failures  of  the  n  units 

on  test,  during  the  time  interval  (0,t].  It  is  assumed  that  {X  (t); 

n 

0  <  t}  is  a  homogeneous  Poisson  process  with  mean  of  X  =  n/6  failures 
per  time  unit,  where  6  is  the  MTBF  of  a  single  unit.  The  graphs  of 
X^(t)  are  non-decreasing  step  functions,  with  jumps  of  size  1  at  the 

random  failure  times. 

The  boundary  line  Lj^(t)  assumes  integer  values  on  [0,t*]  every  A 
units  of  time,  where  A  =  1/b.  Accordingly,  let  t^  =  iA  (i  =  0,l,2,...,kp  + 

The  graph  of  X  (t)  can  cross  the  upper  boundary  b„(t)  at  any  time  t,  but 

X^(t)  can  cross  the  lower  boundary  b^(t)  only  at  the  time  epochs 

t.,  i  >  k^.  Consider  the  time  Interval  (0,t,].  X  (t)  could  have  crossed 

upper  boundary  bj^(t)  at  some  point  1  ,  0  <  T  <  t  <  if  and  only  if 

X^(t)  >  k^  +  1 .  Accordingly,  if  we  denote  by  Pos(j;X)  the  CDF  of  the 

Poisson  distribution  with  mean  X ,  at  a  point  j,  and  by  p(j;X)  the 
corresponding  probability  distribution  function  (PDF);  i.e.,  p(j;X)  = 

-X  i  J 

e  X'^/jl,  end  Pos(j;X)  =  Z  p(i;X),  then 

i=0 

Fg(t)  =  F^^^(t)  =  1  -  Pcs(k^;^),  0  <  t  <  t^.  (5.1) 

Define  the  defective  probability  distribution  function 

g^^(j;t)  =  P0^X^(t)  =  j,  T  >  t).  (5.2) 

-  10  - 
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h-v- 
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It  follows  that,  for  every  j  =  0,...,kj^, 


j  =  0 . k^. 


(5.3) 


Furthermore, 


r-.N-v 
.V.’ 


Pe{T  >  tj)  .  r  8g(tit  )  .  p„s(k  ;  f  t  ).  (5.4) 

To  derive  the  distribution  for  t  values  r.reater  than  t^  we  utilize  the 

property  that  the  Poisson  process  is  an  Increasing  process  of  independent 

2  1  i  1  kjj,  we  apply  the  recursive  equation 


8e(J;t^)  =  J 


•J 

L^o  f  0<j<k^+i-2 

kj^+i-2 

Jlio  j-k^+i-1 


(5.5) 


I 

vK  - 


/-V 

p'-‘ 

PS 


ifP-S 

'mi 

r-T* 


;.v. 


Like  in  (5.4),  for  every  i  *  l,...,kQ, 


PftfT  >  tj  =  Z 


k^+i-1 


i'  - 

j=o 


Furthermore,  for  every  2  <  i  <  k^. 


and 


V'i-1  ‘  i  -  rg{r  >  t._j)  -  P^{7  >  t^) 


k^-hi-2 


(5.6) 


(5.7) 


j^O  f8g(j;ti_i)  -  gQ(j;t.)]  -  g.(k^  +  i  -  l;t^), 

F0(t^)  =  Fe^*^i-i>  +  ^e^*^i-i  ^  1  1  1  i  1  kp  (5.8) 
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where  S  0  and  (0) 
Fg(t^).  However,  at  i 


0.  Notice  that  for  every  i  < 
^0’ 


1,  Fg-'Ct.) 


o  Kq  ®  *^o 


(5.9) 


and 


(5.10) 


t,  is  the  first  jump  point  of  (t).  The  value  of  F^Ct)  at  any 
kQ  be 

t^_j^  <  t  <t^,  i  -  l,...,kQ,  can  be  determined  from  the  previous  formula 
by  substituting  t  for  t^  in  (5.5)  -  (5.8)  and  replacing  A  with 
in  (5.5).  The  corresponding  PDF  of  ^^(t)  can  be  der.Lved  from  the  above 
foimula  by  differentiation.  Alternatively,  we  can  write  for  j 

k^+i"2 


(5.11) 


Pos(k^  +  1  -  1  -  Jl;  I  (t  -  t^_^))] 


Differentiation  with  respect  to  t  yields  the  PDF,  for  ^  *■!’ 

^1^ 


X  k,+i-2 


fe(t)  =  ^ 


e  *^i-i^’^^*^i  ■*■  i  -  1  e(f-  - 


n  ,,  n  . 

e  p<ki:  F  t). 


if  i  =  2  , .  .  .  ,kQ . 


if  i  =  1. 


(5.12) 


We  assumed  that  k  >  k_  +  k, .  Consider  now  the  time  interval 

-  *■  -  *'kg-kj^.  In  this  interval  the  two  boundary  lines  f>j^(i) 
are  parallel.  On  this  interval  we  define  for  every  l<i<k  -k^-k.,, 

*  ^  S  w  X 


12 


T-506 


w’V 


y  1  »  . . 


j 


I  gU;t.  .  ,)p(.i 
£  =  1  0 


i;  A). 


l<j<ko+k^+i"2 


kg+kj^+i-2 


I  g(«.;  t  ^ 

ji=i  V  ^  ^ 


-  1  - 


t;  I  A),  J=kQ+k^+i-l 


The  probability  of  acceptance  at  t  .  is 

kQ+J 


and 


As  before,  we  can  write. 


kp+kj^+i-l 


''e'"  "  Vi’  ■  ®‘^‘‘k„+i’' 

0  j=i+l  0 


and  hence 


kp+k^+i-2 


j=i 


ko+ki+1-1 


j=i+l  ^  V 


^e^^’’^kg+i-l^ 


kQ+ki+1-2 


+  ^  +i-l^  ~  ^e^^’*'k  +i^^ 

j-i+1  *^0  0 


gg(ko  +  k^  +  i  -  i;tj^^^.). 
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But  since  Pg{T^  - 


<  t,  ^.} 

o  Uq+I-I  R  ^0^^ 


(5.18) 


kQ+k^+i-2 


^  l80(jJtj^  +1-1^  ”  +1^^ 

j=i  *^0  0 


-ge(ko  +  k^  +  i  -  i;tj^^^.i). 


Adding  this  increment  to  the  value  of  Fg  (t)  at  t^^  +i_i  obtain  the 
value  of  Fg  (t)  at  tj^  We  can  also  write,  for  every  tj^  +i-i  ^  ^  +i 


Fg’'!!)  .  + 


(5.19) 


ko‘<ki+i-2 


-  "‘’“<‘‘0  +  kj  +  1  -  1  -  f  (t  -  Wi))l 


This  fonnula  yields  by  proper  differentiation  the  PDF,  for 


\*i-i "  ' "  \*i- 


kg+kj+1-2 


ffl(t)  =  ^ 


I  ^e^’^^k.+i-l^ 
£‘=i  0 


(5.20) 


.p(kg-.k^^i-l-£;-(t  -  t,^,,.,)). 


Finally  we  consider  the  Interval  ^  ^  ^k  +k  1’  this  interval 

*^6  *^1  ~  *^s  ^0 


the  upper  boundary  bj^(t)  is  the  constant  k^.  The  functions  gg(j;tj^  ^^) 

6  1 


are  determined  recursively,  according  to  the  formula,  for  every 


^s  -  >^1  -  *^0  ^  i  ^  J  ^  ^ 


-  14  - 


-  V 


I 

4 
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-k,+i>  ■  .  ^  ^  8e“’"k  -k.+i-i’fO  - 

s  1  J!,«=k  -k,-k-+l  si 

s  1  u 


l<l<ko+ki-l. 


From  this  one  obtains  immediately  that 


"  \  -k  +i^  "  ■  *^0  ■  ^1  ■*■  ^’‘^k  -k  +i^’ 

si  si 


^e^'^k  -k  +i-l  ^  ^  -  ‘^k  -k  +i^ 

si  si 


^3^*^s  ■  ^0  "  **^1  ^’*^k  -k  +i-l^ 

6  1 


s  0  1 


Accordingly, 


^*^k  -k  +i^  ^  ^*^k  -k  +i-l^  ''■ 

si  si 


^e^'^s  “  *^1  “  ^0  '*’  ^’*^k  -k,+l^’ 

s  1 


^‘^k  -k  +1^  "  ^*^k  -k  +i-l^  ■'■ 

si  si 


(5.21) 


(5.22) 


(5.23) 


(5.24) 


(5.25) 


i!^k  -k  -k.+i 
E  1  0 


-k  +i-l^  ■  -k  +i^^* 

si  si 


-  15  - 


■■  v' 


-  ■  -  •  ».iVj  W»V,.  aV  it\.  «*A  »  _  if.  »  j  J".  f.  .' 


T-506 


For  values  of  t  in  the  Interval  (tj^  +i-l‘*^k  -k  obtain 

*^s  1  ^1  ^ 


S  1 


^  ,)(1  -  Pos(k  -  1  -  £;  ~(t  -  t,  ,  , 

£=k6-ki-ko+i  ®  ®  ® 


fe(t) 


ks"! 

t  ^  -k  +i-l^' 

t-k^-kj-kfl+l  '‘b  '‘1+^ 


•p(k^  -  1  -  t:  e  (t  -  tk  .. 

S  1 


6.  The  Expected  Stopping  Time  (ATL) 

The  expected  stopping  time,  or  average  test  length  (ATL)  is 
given  by 


^  t  P  {t  =  t  ] 

1=0  ®  A 


‘'o'"‘‘b-*  fK. 

"  ) 

V=1  •'t 


tdF  ''(t). 


The  first  term  on  the  r.h.s.  of  (6.1)  is  equal  to 


'k.+i  8e“-‘k>i> 
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The  second  term  on  the  r.h.s.  of  (6.1)  is  expressed  as  a  sum  of  three 
R  R  R 

components,  E,  q  »  E.,  ^  and  £_  _  ,  where 

X|U  ^|U  «5^t7 


0  A 

I  f 


1  K-i 


tdFg"(t). 


k  -k,  t 

_  6  1  r  V  _ 

E;  -  I  tdY^it), 

"■®  -VI  ’vi  ' 


(6.3) 


k  +k.-l  t 

s  0  r 


V=k  -k,+l 
s  1 


\  tdF  ^(t). 

•^Vi  ® 


(6.4) 


The  following  result  is  applied  in  the  development  of  the  formulae  for 
E^^d  =  1,3);  namely. 


I  tp(j;X(t  -  t^_^))dt  = 

Vl 

=  -Y^  0  +  1)(1  -  Pos(j;XA))  4  (1  -  Po3(j  +  1;XA)). 


(6.5) 


According  to  (5.12)  and  (6.5), 
ko  k^+1-2 

•(1  -  Pos(kT  +  i  -  £  -  1;  §  A))  +  ^  (k,  +  i  -  £)(kT  +  i  +  1 
i  W  n  J.  i 


(6.6) 


(1  -  Pos(k^  +  i  -  £;  I  ))}. 
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(t 


mourn 


For  Q  we  use  (5.20)  and  (6.5)  and  obtain  the  formula: 


k  -k  -k^  k-+k,+i-2 
s  1  0  0  1 


{tk^^.i_l(ko  +  k^  +  1  -  £)[1  -  Pos(kQ  +  k^  +  1  -  1  -  I  A)]  (6.7) 


+  ~  (k^  +  kj^  +  i  -  2^)(kQ  +  kj^+±  +  l-  £,)[!  -  Pos(kQ  +  k^^  +  i  -  £  A)  ] } 


Finally,  from  (5.27)  and  (6.5)  we  obtain 


k-+k,-l  k  -1 

0  1  8 


i-1  £.=k  -k  -k.+i 

6  1  0 


-k,+i-i^ 

s  1 


^‘■k  -k,+i-i^^s  ■  ^ 

s  1 


f  •  (k^  -  £)(k^  +  1  -  £)[1  -  Pos(k^  -  I  A)]} 


Formulae  (6.2),  (6.6).  (6.7)  and  (6.8)  yield  the  expected  stopping  time. 


7.  Concluding  Remarks 


As  shown  in  the  previous  sections,  in  order  to  determine  a  lower 


confidence  limit  for  the  MTBF,  6,  after  stopping  according  to  the  rule.s 


of  sequential  life  testing,  one  has  to  graph  the  percentile  points 


Ta  distribution  F^  (t)  as  a  function  of  6.  This  requires 


A  R 

the  tabulation  of  the  distributions  F„  (t),  F„  (t)  and  FQ(t),  for  various 

DO  D 


values  of  0.  In  the  Appendix  we  provide  a  BASIC  program  for  the  computa- 


R  A 

tion  of  the  cumulative  distributions  F„  (t),  F^  (t)  and  FQ(t)  for 

U  O  V 


various  values  of  0.  After  compilation,  it  has  taken  about  five  minutes 
to  run  this  program  on  an  IBM  PC,  with  five  different  values  of  6. 
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By  changing  the  parameters  k_,  k  ,  k  and  b  of  the  stopping  boundaries 

O  X  S 

one  changes  the  characteristics  of  the  procedure.  Thus,  if  it  is  desired 
to  change  the  acceptance  probabilities  tt(6)  or  the  average  test  length, 
ATL(6),  it  is  possible  to  try  different  values  of  the  boundary  line 
parameters  and  examine  the  resulting  graphs  of  7r(0)  and  ATL(0). 


f  V'.3 


'J  •W^TTTrW 
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APPENDIX 


BASIC  PROGRAM  FOR  THE  COMPUTATION 
OF  THE  DISTRIBUTIONS  OF  STOPPING  TIMES 


10 

PRT 

NT  " 

A 

PROGRAM 

20 

PRT 

NT  " 

i 

nsert 

va 

30 

INP 

UT  N,k 

0 

,K1 

,KS 

!.E< 

40 

KB= 

KS-KO- 

K 

1 

50 

KBP 

=KB+1 : 

K 

BM= 

KB- 

1 

60 

l:  1  p 

=K1  +  1 : 

k: 

1M= 

Kl- 

1 

65 

l::sp 

=KS-t-l : 

k: 

SM= 

k:s- 

1 

70 

DLT 

=  1/B 

75 

DIM 

H(KSF 

) 

,G( 

KSP 

)  , 

FOR  COMPUTING  THE  CHARACTERISTICS  OF  THE 
ues  o-f  n,  kO,  PI ,  I:s,  b  in  order" 


PC (KS+KO) 


80  TETJ=3500 

90  FOR  TET=1500  TO  5500  STEP  1000 
100  W=TET/TET1 


DESIGN 


no  TO=KOU(DLT 

120  LPRINT  USING"  ##  ##  ##  ##.♦♦  ;  N,  KO,  K1  ,  KS,  B,  TET 

130  HET=DLT»N/W 

135  LAM=HET 

140  FOR  1=0  TO  KS 

150  IF  I>  K1  THEN  GOTO  230 

170  M=1 


180  GOSUB  1300 
185  FI=PS 
190  M=M-1 
200  GOSUB  1300 
210  H(I)=FI-PS 
220  GOTO  240 
230  H(I)=0 
240  NEXT  I 


250  M=K1 
260  GOSUB  1300 
270  CDF1=1-PS 
280  PC ( 1 ) =CDF 1 
290  Tl=Di-T 

300  LPRINT  USING  "  #.#«##";  T1,CDF1 

310  FOR  1=2  TO  kO 

320  ^..I=^  1  +  I 

330  k;ip=k:i  +  i :  l:im=li-i 

340  FOR  J=0  TO  LIM 

350  SUMJ=0 

360  FOR  L=0  TO  J 


370  JL=J-L 

380  GOSUB  1300 

390  PI=PS 


400  M=M-1:  GOSUB  1300 
410  SUMJ=SUMJ+H (L) * (FT -PS) 


420  NEXT  L 
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450  G(J)=SUMJ 
460  NEXT  J 
4X0  FOR  J=KI  TO  KS 
480  G(J)=0:  NEXT  J 
490  TI  =  I)t:DLT 
500  TUM=0 

510  FOR  L=0  TO  (KIM-1) 

520  TUM=TUM+H(L)-G(L) 

530  NEXT  L 

540  CDF1=CDF1+TUM-G (KIM) 

550  PC(I)=CDF1 

555  IF  I=KO  THEN  GOTO  570 

560  LPRINT  USING"  I).####  •■;TI,CDFl 

570  FOR  J=0  TO  KS 
580  H(J)=G(J) 

59C>  NEXT  J:NEXT  I 

600  PT(0)=H(0) :CDF2=PT(0) 

610  CDFT=CDF1+CDF2 

620  LPRINT  USING  "  #.♦»#  #.###«  ♦I.##**#  T  I ,  CDF  1  .  CDF2 ,  CDFT 

630  FOR  1=1  TO  KB 

640  FOR  J=0  TO  I 

650  IF  J<I  THEN  6(J)=0 

660  NEXT  J 

670  KI=K0+K1+I:  KIM=KI-1 
600  FOR  J=I  TO  KIM 
690  SUMJ=0 
700  FOR  L=I  TO  J 
710  JL=J-L 

720  M=JL:  60SU&  1300 
730  PI=PS 

740  M=M-1  :GOSUB  1300 
750  SUMJ=5UMJ+H (L)  *  (PI-PS) 

760  NEXT  L 
770  G(J)=SUMJ 
775  NEXT  J 
780  PT ( I ) =G  ( I ) 

790  CDF2=CDF2+PT (  1  ) 

80<.>  TI=T1+DLT 
810  TUM=0 

820  FOR  L  =  I  TO  ()  IM-1  ) 

B3<:>  IF  L=I  THEN  TUM  =  TUM  +  H(L)  ELSE  TUM=TUM  +  H  ( L  ) -G  ( L  ) 

840  NEXT  L 

850  CDF1=CDF  1+TUM-G  O  IM) -PT  (  I  ) 

860  PC (KO-H  ) =CDF1 :  CDFT=CDF 1 +CDF2 

87'  ’  LPRINT  USING  '  T  1  .  CDF  1  ,  CDF2,  CDFT 

880  FOR  J=0  TO  F5 

890  H(J)=G(J):  NEXT  JrNEXT  I 

900  FOR  I=KBP  TO  fSM 


22 


910  FOR  J=0  TO  (I-l) 

920  G(J)=0  :NEXT  J 
9iO  FOR  J=I  TO  KSM 
9 AO  BUMJ=0 
950  FOR  L=I  TO  J 
960  JL=J-L 

970  M=JL:  GOSUB  1300 

980  F'I=PS;  M=M-1:  GOSUB  1300 

990  SUMJ=SUMJ+H (L) * (PI-PS) 

1000  NEXT  L 

1010  6(J)=SUMJ  :NEXT  J 
1020  PT(I)=G(I) 

1 030  CDF2=CDF2+PT ( I ) 

1040  TI=TI+DLT 

1050  Turi=o 

1060  FOR  L=I  TO  KSM 

1070  IF  L=I  THEN  TUM=TUM+H(L)  ELSE  TUM=TUM+H (L) -G (L) 

1080  NEXT  L 

1090  CDF1=CDF1+TUM-PT(I) 

1100  PC  (KO+I )  =--CDFl :  CDFT=CDF1 +CDF2 

1110  LPRINT  USING  "  #.##  #.####  #.####  #.####"; T I . CDF  1 , CDF2 . CDFT 

1120  FOR  J=0  TO  KS 

1130  H(J)=G(J) :NEXT  JrNEXT  I 

1140  ASN=0 

1150  KSV=KS+KO 

1160  FOR  1=1  TO  (KSV-1) 

1170  TI=I#DLT:  LI=1-K0 

1180  IF  I>=  KO  THEN  ASN=ASN+TIirPT  (LI ) 

1190  IF  1=1  THEN  PF=PC(1)  ELSE  PF=PC ( I ) -PC (I - 1 ) 

1200  ABN=ASN+ (TI-DLT/2) *PF 
1210  NEXT  I 
1220  DCP=CDF2 

1230  LPRINT  USING  "  ocp=#.##tttt  a5n=« . ; OCR . ABN 

1235  NEXT  TET 
1250  END 

1300  IF  M<0  THEN  PS=0 

1310  IF  MCO  THEN  GOTO  1400 

1320  PR=EXP (-LAM) 

1330  IF  M=0  THEN  PS=PR 

1350  IF  M=0  THEN  GOTO  1400 

1355  CP=PR 

1360  FOR  KP=1  TO  M 

1370  PR=PR)rLAM/KP 

1380  CP=CP+PR 

1390  NEXT  KP 

1395  PS=CP 

1400  RETURN 


OFFICE  OF  NAVAL  RESEARCH 
STATISTICS  AMD  PROBABILITY  PROGRAM 

BASIC  DISTRIBUTION  LIST 
FOR 

UNCLASSIFIED  TECHNICAL  REPORTS 


FEBRUARY  1932 


Copies 


Statistics  and  Probability 
Program  (Code  411 (SP)) 

Office  of  Naval  Research 
Arlington.  VA  22217  3 

Defense  Technical  Infcrmatlon 
Center 

Cameron  Station 

Alexandria.  VA  22314  12 

Commanding  Officer 
Office  of  Naval  Research 
Eastern/Central  Regional  Office 
Attn:  Director  for  Science 
Barnes  Building 
495  Summer  Street 
Boston,  MA  02210  1 

Commanding  Officer 
Office  of  Naval  Research 
Western  Regional  Office 
Attn:  Dr.  Richard  Lau 
1030  East  Green  Street 
Pasadena,  CA  91101  1 

U.  S.  ONR  Liaison  Office  -  Far  East 
Attn:  Scientific  Director 
APO  San  Francisco  96503  1 

Applied  Mathematics  Laboratory 
David  Taylor  Naval  Ship  Research 
and  Development  Center 
Attn:  Mr.  6.  H,  Glelssner 
Bethesda,  Maryland  20084  1 

Conmandant  of  the  Marine  Coprs 
(Code  AX) 

Attn:  Dr.  A.  L.  Slafkosky 
Scientific  Advisor 
Washington,  DC  20380  1 


Copies 


Navy  Library 

National  Space  Technology  Laboratory 

Attn:  Navy  Librarian 

Bay  St.  Louis,  MS  39522  1 

U.  S.  Anry  Research  Office 
P.O.  Box  12211 
Attn:  Dr.  J.  Chandra 
Research  Triangle  Park,  NC 
27706  1 

Director 

National  Security  Agency 

Attn:  R51 ,  Dr.  Maar 

Fort 'Meade,  MD  20755  1 

ATAA-SL,  Library 
U.S.  Army  TRADOC  Systems 
Analysis  Activity 
Department  of  the  Army 
White  Sands  Missile  Range,  NM 


88002  1 

ARI  Field  Unlt-USAREUR 

Attn:  Library 

c/o  ODCSPER 

HQ  USAEREUR  &  7th  Anny 

APO  New  York  09403  1 

Library,  Code  1424 
Naval  Postgraduate  School 
Monterey,  CA  93940  1 

Technical  Information  Division 
Naval  Research  Laboratory 
Washington,  DC  20375  1 

OASD  (I&L),  Pentagon 

Attn:  Hr.  Charles  S.  Smith 

Washington,  DC  20301  1 


Copies 


Copie 


D1 rector 
AMSAA 

Attn:  DRXSY-MP,  H.  Cohen 
Aberdeen  Proving  Ground,  MD  1 
21005 

Dr.  Gerhard  Helche 
Naval  Air  Systems  Cotnnand 
(NAIR  03) 

Jefferson  Plaza  No.  I 
Arlington,  VA  20360  1 

Dr.  Barbara  Ballar 
Associate  Director,  Statistical 
Standards 
Bureau  of  Census 
Washington,  DC  20233  1 

Leon  Slavin 

Naval  Sea  Systems  Command 
(NSEA  05H) 

Crystal  Mall  #4.  Rm.  129 
Washington,  DC  20036  1 

B.  f.  Clark 

RR  #2,  Box  647-8 

Graham,  NC  272S3  1 

Naval  Underwater  Systems  Center 
Attn:  Dr.  Derrlll  J.  Bordelon 
Code  601 

Newport,  Rhode  Island  02840  1 
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